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We investigate the stability of large-scale structures in Burgers' equation under the 
perturbation of high wave-number noise in the initial conditions. Analytical estimates 
are obtained for random initial data with spatial spectral density < 1. Numerical 
investigations are performed for the case n = 0, using a parallel implementation of the 
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The appearance of ordered structures is possible for nonlinear media with dissipation 0] . 
These structures often form a set of cells with regular behaviour, alternating with ran- 
domly localised zones of dissipation. One nonlinear dissipative system with such behaviour 
is the well-known Burgers' equation 

dv dv d'^v f \ f \ /X 

Tt^'Tx=^^-'^ .(x,t = 0)=.o(x) (1) 

where /i is the viscosity coefficient. The solutions of Burgers' equation with random initial 
conditions display two mechanisms, which are also inherent to real turbulence: nonlinear 
transfer of energy through the spectrum, and viscous damping in the small scale region. 
Equation (0) was proposed by Burgers^] as one-dimensional model of fluid turbulence. It 
has since been shown to arise in a large variety of non-equilibrium phenomena, when parity 



invariance holds|TO|. Burgers' equation have applications to nonlinear acoustics, nonlinear 
waves in thermoelastic media, to modelling of the formation of large-scale structures in 
the Universe , and to many other systems where dispersion is negligibly small compared 
with nonlinearity (see and references therein). 

Burgers' turbulence is an example of strong turbulence, the properties of which are 
determined by the strong interaction of large number of harmonic waves. In Burgers' 
turbulence, due to the coherent interaction of harmonics, saw-tooth shock waves are 
formed, which may be treated as a gas of strong local interaction of particles [l^, H . The 



collisions of shock fronts leads to their merging, which is analogous to inelastic collisions 
of particles, and to an increase of the integral scale of turbulence. 

In the limit of zero viscosity the solution of (|I]) for the velocity field is given by 

x-y{x,t) dS{x,t) 
v{x,t) = ; v{x,t) = — — (2) 

where y{x,t) is where is realised the absolute maximum of the function 

G{x,y,t) = So{y) — — ; So{x) = J vo{y)dy (3) 

By analogy between the solutions of Burgers' equation and the fiow of particles P , we 
shall call So{y) the initial action, and S{x,t) = G{x,y{x,t),t) the action at time t. This 
solution can be computed by a Legendre transform of initial data |T2[. We shall assume 



that the spectrum of random initial velocity field has a following form: 

goik) = alk^h^{k) boik < k^) = 1, ho{k > fco) = (4) 
The curvature Kg of the initial action in may then be estimated as 

Kl =< {SMr >= c^lkrVin + 3) = aoV/Ln, =< >= "n^o^Vl^ + 1) (5) 

Taking into account that the curvature of parabola in is 1/t, we get that when 
t 3> = ^mm/co = the curvature of the parabola is much smaller than the curvature 
of the initial action. This implies that the global maximum of G{x, y, t) is in the neigh- 
bourhood of the local maximum of So{y), and that y{x,t) is a stepwise non-decreasing 
function of x. Thus the velocity field has a universal behaviour 

v[x,t) = — - — ; Xk^i<x<Xk (6) 
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in each cell between shock positions Xk- Each cell is characterized by two numbers: the 
inverse Lagrangian coordinate and the value of initial action in the cell Sk = So{yk)- 
The dissipation zone is in this case at the shock front, and hence have a zero width. The 
positions of the shocks are determined by the equality at the point Xk of the values of two 
absolute maxima G{xk,yk-i,t) = G{xk,yk,t): 

Vk + Vk-i , ... So{yk-i) - So{yk) , . 

Xk = h Vkt; Vk = (7) 

2 yk- yk-i 

It is easy to see that the rate of collisions of shocks depends on the asymptotic behaviour 
of the structure function of initial action 



4(p) =< [Soix + p) - iSoix)Y > 



alp^ if < 1 
2al if n > 1 



Here =< Sq >. We can estimate the integral scale of turbulence l{t) = \x — y\ from 
the condition that the parabola and the initial action are of the same order: 

dsimY^'^iVt (9) 

From it therefore follows that we have two different types of growth for the scale 

i{ty. 

r («t)V("+3)^ if-<i .,0. 



It was shown |^ (see also [|T2[) that the common feature for both types of initial spectra, is 
the existence of self-similarity of statistical properties of solutions, which are determined 
by only one scale l{t). From (|TUp one can see that in the case n < 1, the behaviour of 
Burgers' turbulence at time t is determined by the local large-scale behaviour of the initial 
spectrum. This was the motivation for a simple qualitative model of Burgers' turbulence 
(at n < 1) as a discrete infinite set of modes - the spatial harmonics km = fco7~'"(7 ^ 1)? 
sufficiently spaced in the spatial spectrum ^ ||. The amplitudes Af^ of harmonics were 
choosen from the condition that the mean spectral density of harmonics in the interval 
= km+1 — km was identical to the random noise spectral density: = gQ{km)^m- 
Assuming that the influence of large-wavenumber components on large-scale one is small, 
it was obtained that the growth of average scale for this discrete model is the same as 
in the case of continuous spectrum given by (|1^). The idea of independent evolution of 
large-scale structures has also been used to obtain long time asymptotics of the cylindrical 
Burgers' equation [Q. 

It is possible to show that the assumption of relatively independent evolution of the 
large-scale structures is valid for the continuous model of Burgers' turbulence as well. 
In it was obtained analytically that the interaction of a regular positive pulse with 
noise in Burgers' equation does not change the evolution of the large-scale structure of 
the pulse, when the value of initial action of the pulse is greater than the dispersion of the 
noise action. This result was confirmed by numerical simulations based on the asymptotic 
solution (^H) of Burgers' equation. Here we shall give some simple estimations of this 
effect for the random perturbations with initial spectrum (^) and illustrate it by results 
of numerical experiment for the case of truncated white noise (n=0 in (pUj)). 

Numerical experiments were performed using a parallel version of Fast Legendre Trans- 
form alghorithm implemented on a Connection Machine CM-200 H 0. The algorithm 
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uses the property that the inverse Lagrangian function y{x, t) is a non- decreasing function 
of X, and specific low-level instructions on the Connection Machine, that allows one to 
perform simultaneously partial maximization operations over divisions of the range, that 
are only known at run-time. The most primitive computation of the maximum of (^) 
takes for points 0{N'^) operations on a sequential computer. The serial Fast Legendre 
Transform takes 0{N\ogN) operations in ID (and 0{N'^ log^ A^) operations for data on a 
N X N grid in 2D [|11|]). Our Parallel Fast Legendre Transform takes 0(log A^) operations 
on an ideal parallel computer with unlimited number of physical processors, connected in 
hypercube network (this is the type of connection used on the CM-2 and CM-200 mod- 
els). The Connection Machine simulates an arbitrary number of processors with a finite 
number (on our machine 8192). Our algorithm will therefore eventually go linearly 

with the ratio of initialized gridpoints (A^), to the actual number of physical processors. 
For the class of one- dimensional problems we have considered here, the computational 
time including input, output and Connection Machine initialization is much less then one 
minute for total number of points up to A^ = 2^°. 

We shall consider the evolution of two initial random perturbations Vo{x) and ^0(3;): 

v^{x) = Vo{x) + Vhix) (11) 

The power spectra of both processes are described by (^ with n < 1, but the process 
Vo{x) has spectral components in the range [0,A;*], while the process Vq{x) in the range 
[0, ko] with fco ^ fc*. Restricted to the range [0, k^] the two processes are identical, Vh{x) 
being their difference in the large wave-number range [k^, ko]. The correlation coefficient 
between vq and Vq is tq = {k^/koY'^^^^^'^ <^ 1. In particular, we made calculations for 
JvQ ^ ^ ^ 2 cLiid A^Q ^ I k-^ ^ 2^. 

Three initial realisations, as described above, are shown on Fig.l. At t ^ t„ = 
l/(aA;i"+')/') all three processes are transformed to a sequence of triangular pulses with 
the universal behaviour inside the cells according to (|^). The solution f(x,t) will be 
stable relative to high wave- number perturbation Vh{x), if the fluctuations of both the 
inverse Lagrangian coordinates Ayt = yk — yk and the shock positions Axk = Xk — Xk, are 
small respectively to integral scale of turbulence l{t). While the asymptotic properties of 
Burgers' turbulence are determined by the behaviour of initial action, the values of dis- 
turbances Axk and Ayk will be determined by the value of dispersion a^g of perturbation 
action Shix) = ^ Vhiy)dy: 



a? 



l-nN 



'^-=<^'>=(r3^|i-|T^l I (12) 

and its correlation scale Ig ~ (/c^,/co)^^^^. The main divergence between v and v is due to 
the different velocities of shocks Vk and Vk, and these errors increase in time. From (|^) 
one sees that the fluctuations of velocity are determined by the fluctuations of the action 
Sh{y) in the neighbourhood of local maximum of Soly). This problem is similar to the 
analysis of statistical properties of Burgers' turbulence with stationary initial action 
Thus we get that the relative deviation of shock positions is equal to 

eit) - < ^ ^ (13) 
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which means that the behaviour of e{t) depends on the growth law of the integral scale 
l{t). From (|rU|,|T2|,|T^ we get that in the case n < 1 and k^: 



Here we have taken into account that for k^ ^ A;*, the value of absolute maximum 
of the random process Sh{y) in the vicinity of local maximum of So{y), has a double 
logarithmic distribution with average value ~ ahs{^T^{ko/ K))^^"^ , and dispersion decreasing 
with correlation scale Is, proportional to ln(A;o/A;*) (|^ ^^)- K may be that one shock in 
V corresponds to a cluster of shocks in v, but then they are spread out over a distance 
no more than e (IT^)- Therefore one can say that the large scale structures of Burgers' 



turbulence are self-preserving due to multiple merging of shocks. The stability of large- 
scale structures is illustrated by results of numerical calculations, presented on Fig. 2. 
Three realisations of velocity field are shown, corresponding to ininitial conditions on 
Fig.l with cutoff wavenumbers equal respectively A;*, k^^^ and kl^\ One can see that the 
large-scale behaviour of all those realizations is similar, and only fine structure weekly 
depends on cutoff wavenumber. The correlation coefficient between the velocity fields 
with cutoff wavenumbers k^ and k^^^ increase from the value 0.035 for initial perturbation 
(Fig. 1) up to the value 0.92 on the stage of developed shocks (Fig. 2). This become even 
more clear when comparing the corresponding initial and final actions, presented on Fig. 
3 and Fig.4. There one sees that the action for higher cutoff values can be obtained from 
the smaller one approximately by a shift along the vertical axis. This is in accordance 
with estimations following from double logarithmic distribution of Sh{y)- 

In conclusion, we would like to stress that the stability of large-scale structures rel- 
ative to large wave-number perturbations, is similar to the independence of asymptotic 
evolution of Burgers' turbulence (with n < 1) of the value of viscosity coefficient /ip. 
We may thus interprete the action of small-scale perturbations on solution of Burgers' 
equation as a turbulent viscosity, effectively changing /i, but leaving the initial conditions 
invariant M . 
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Fig. 3. Realizations of the action fields with cutoff 

wave-numbers k^{l) , kQ^\2) and k^\3) at 
t^tn and initial action for wave-number 
k = fc*(4) . 
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Fig. 4- Realizations of the action fields with cutoff 

(2) 

wave-numbers k^{l) and (2) att^t^ 
and initial actions for the same wave-numbers 

(3,4). 



